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Abstract 



from the nucleoli matrix element of the hyperon correlation function. Through 



J> ' We investigate the hyperon-nucleon interactions in the QCD sum rule starting 

cn ■ 
(N 
O 
00 

(•~^ ' the dispersion relation, the correlation function in the operator product ex- 

cn ' 

' pansion (OPE) is related with its integral over the physical energy region. 

The dispersion integral around the hyperon-nucleon (YN) threshold is iden- 
^ ' tified as a measure of the interaction strength in the YN channel. The Wilson 



coefficients of the OPE for the hyperon correlation function are calculated. 



^ ' The obtained sum rules relate YN interaction strengths to the nucleon ma- 



trix elements of the quark- gluon composite operators, which include strange 
quark operators as well as up and down quark operators. It is found that 
the YN interaction strengths are smaller than the A^A^ interaction strength 
since the nucleon matrix elements of strange quark operators are smaller than 
those of up and down quark operators. Among YN channels AA^ channel has 
stronger interaction than T,N and 3N channels. Also found is that the in- 
teraction strength is greater in the Ti^N (B^N) channel than in the Ti^N 
(H~A^) channel since the nucleon matrix elements of up quark operators are 
greater than those of down quark operators. The spin-dependent part is much 
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smaller than the spin-independent part in the YN and NN channels. The 
results of the sum rules are compared with those of the phenomenological 
meson-exchange models. 
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I. INTRODUCTION 



One of the most important goals of studying strong interaction physics is to understand 
the behaviour of hadrons and hadronic interactions on the basis of the quantum-chromo- 
dynamics (QCD). One of the powerful tools for this purpose is the QCD sum rule proposed 
by Shifman, Vainshtein and Zakharov [1]. Using the method of the QCD sum rule one can 
obtain relations between the properties of hadrons and the vacuum matrix elements of the 
quark-gluon composite operators [1,2]. 

The idea of the QCD sum rule has been extended to the investigation of hadronic inter- 
actions [3-7]. In particular, in Ref. [7] the present authors have formulated the QCD sum 
rules for spin-dependent nucleon-nucleon interactions and have studied their physical impli- 
cations. The basic object of the study in the work has been the nucleon matrix element of 
the correlation function of the nucleon interpolating field. In the deep Euclidean region the 
correlation function is expressed in terms of the nucleon matrix elements of the quark-gluon 
composite operators using the operator product expansion (OPE), which is related with its 
integral over the physical region by means of the dispersion relation. The dispersion integral 
of the correlation function around the nucleon threshold has been investigated in detail. 
It has turned out that the integral can be identified as a measure of the nucleon-nucleon 
interaction strength, which is proportional to the scattering length in the small scattering 
length limit and to the one half of the effective range in the large scattering length limit. 
The obtained sum rules relate the nucleon-nucleon interaction strengths with the nucleon 
matrix elements of the quark-gluon composite operators. The sum rules tell us that the 
interaction strength in the spin-singlet channel is weaker than in the spin-triplet channel, 
but that the spin-dependent part of the interaction strength is considerably smaller than the 
spin-independent part. Experimentally, it has been found that there is a loosely bound state, 
deuteron, in the spin-triplet channel, while there is an almost bound state in the spin-singlet 
channel, which implies that the interaction is slightly stronger in the spin-triplet channel 
than in the spin-singlet channel. This is consistent with the result of sum rules. 
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It is straightforward to apply the method of Ref. [7] to hyperon-nucleon channels just 
by replacing the nucleon interpolating field by the hyperon interpolating field. The purpose 
of the present paper is to carry out this procedure. Physically, it is extremely important 
to understand the various hyperon-nucleon and nucleon-nucleon interactions simultaneously 
in order to achieve a unified view of hadronic interactions. It is also important in the 
application to hypernuclear physics since hyperon-nucleon interactions are exprimentally 
not well known. The paper is organized as follows. In Sec. II, we review the formalism in 
the hyperon-nucleon (YN) channel for the sake of completeness. In Sec. Ill, we calculate 
the Wilson coefficients of the OPE for the hyperon correlation function and explain how 
the nucleon matrix elements of the quark-gluon composite operators are determined. The 
obtained sum rules for the hyperon-nucleon interaction strengths are presented and discussed 
in Sec. IV. Finally, we summarize the paper in Sec. V. 

II. PHISICAL CONTENTS OF CORRELATION FUNCTIONS 

We first consider the spin-dependent hyperon correlation function, Il{q'ps), 



where \ps) is the one-nucleon state with momentum p and spin s {p^ — Mff, — —1 and 
ps — 0, where M^v is the nucleon mass) normahzed as {ps\p's') — (27r)^5^(p — p')Sss' and ■0 
is the normalized hyperon field operator. In this paper, momentum with " represents the 
on-shell momentum. Later, the normalized nucleon field, ijj, is replaced by the unnormal- 
ized hyperon interpolating field (quark-gluon composite field), rj. The following discussion, 
however, holds as it is for the interpolating field, except for the normalization. Naively, the 
dispersion relation for the correlation function, Il{qps), is written as 



where q' = {q^, q). Throughout this paper, whenever we take the imaginary part of a quan- 
tity, we approach the real energy axis from above in the complex energy plane. Therefore, 
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(1) 



strictly speaking, Imll is the imaginary part of the retarded correlation function. The QCD 
sum rules are obtained by evaluating the left-hand side of Eq. (1) by the OPE and expressing 
the right-hand side in terms of physical quantities. 

Let us consider the singularities of Il{qps) as functions of go- In the complex go plane, 
n(gps) has a branch cut from the lowest B — 2 continuum threshold to the right and another 
branch cut starting from the lowest B = continuum threshold to the left where B denotes 



the baryon number, n(gps) has second-order poles at qo = + M"^ = ^Eq, where M 

denotes the hyperon mass. In addition, the coefficient of the pole at go = Eg is the YN 
T- matrices T+: 

Af 

T+{qrps; qrps) = {qo - Eqf —u{qr)-foIl{qps)jou{qr) 

^ qo=Eq 

where u{qr) is a positive energy solution of the free Dirac equation for the hyperon. 

In order to take out the pole contribution from Imn(gps) it is convenient to define off- 
shell YN T-matrices by 

T^{q'r'ps'; qrps) 



^-ij dV-^^u{,ir'){i' - M){p's'\T{^ix)m)\ps){i - M)^^u{qr). (2) 

Note that Eq. (2) is just a definition of the T-matrix off the mass shell, but the LSZ reduction 
formula shows rigorously that it is the T-matrix on the mass shell. 

In order to eliminate the pole contribution at go = —Eg from the correlation function, 
we introduce the projection operators A_|_ by 

f + M 1 + 75/ 



A+(gr) = u{qr)u{qr) = 



which have the properties 



Aliqr) = A+{qr), 
7oA+(gr)7o = A+(gf). 



where g = (go, —q) and r = (— ro, r) 



Then we define the projected correlation functions by 

n+{qrps) = ^ti {A+{qf)U{qps)} . (3) 
The projected correlation functions are related to the off-shell T-matrices as 

where T+{qrps) = T+{qrps; qrps). Clearly, Il+{qrps) has a second-order pole at go = Eg but 
not at qo = —Eg. 

We next consider the dispersion relation for 'n.^{qrps) , 



— —U+{qrps) = -- [ dq'g Im i ^n+(gVps) I . 

qo TT j-oo qo-qo + ^1^ [ qo J 



(4) 



Symmetrizing Eq. (4) we obtain 

— U+{qrps) + {qo -qo) = - dq'^- — — — ^ ^(qq - Eg)lvDli+{q'rps), (5) 

^qo ^ J—oo \qO "r IT]) 

where 

p 

{qo - Eg)lmli+{qrps) = -Tr6{qo - Eg) ReT+{qrps)\ H —lmT+{qrps). (6) 

y qo — 

Applying the Borel transformation, 

Urn ( 



to both sides of Eq. (5), we obtain 
Lb 



^° ^'^U+{qrps) + {qo -qo) 



2qo 



where Mb is the Borel mass. In order to derive the Borel sum rules we must evaluate 
the left-hand side by the OPE and parameterize the right-hand side in terms of physical 
quantities. 
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We study the behaviour of the two terms on the left-hand side in Eq. (6). For this 
purpose it is important to note that the off-shell optical theorem holds for T. When the 
center-of-mass energy is above the threshold of the YN channel and below the threshold of 
the next channel, only the YN states contribute in the intermediate states, and the off-shell 
optical theorem is simplified as 



(27r)3 (27r) 

(8) 



7^^7^^(27r)^(5^(p + q-Pn- qn)T+{qp; qnPn)T+{qnPn] 



In order to simplify the notation we introduce the scattering amplitude / by 

f{q'p';qp) = — T+{q'p';qp), 

where /i = and n' — . Moreover, we go to the center-of-mass frame {q + p = 

q'+p' — 0) and restrict ourselves to the s-wave. We define the off-shell scattering amplitudes 
as 

f(k)^f(q'p';qp), 

where \p\ = \q\ = \p'\ = |qf'| = 0, po = Po = Mjv, ?o = = VM^ + + Ml + k'^ - M, 
From Eq. (8) we have 

Im-^ = -A; + 6A;3 + 0(P), 

fyk) 

and find that the off-shell scattering amplitude, /(/c), has the form 

1 



N- 



fik)^ 



i {-k + bk^ + 0{k^)} + + Ifk-" + 0(A;4)} ' 

where a is the scattering length. It should be noted that f is different from the effective 
range r, but f coincides with r in the limit a — > oo: 

f — r -\- O ^-^ . 

We proceed to the integral of the right-hand side of Eq. (7), /, in the vicinity of go — M, 



The integral, /, can be decomposed as 

I^It + Ic{+h). (9) 
In Eq. (9), the first term, /j, is the threshold contribution, given by 

The second term, I^., is the continuum contribution, given by 

^12 



where 



Im/I"* = ^ 7 ^Ax ^(?o - M). 

^ + + ^ + + 

The last term, J^, is the bound-state contribution, which has to be taken into account if 
there is a bound state, given by 



where 



T rpole 



9 /l^ 



= 7rc5{K — kq), 



S{k- ko) 



k=iKo ■ 



where k — ik and Ikq is the pole momentum, 1/ f2{iK,o) — 0. 

By performing the integral, the continuum contribution becomes 



yjv _ -L I M'^\M + Mn 2n\a\ 



tYN 



Ml MMn /i I f I 6 I 1 , , ~' 

Y a a2 2a2M2 ^ Mat; '^^M^ 

which is simplified in two limits of a as 
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rYN 



Mi 



■ exp 



1 



M2 

Ml 



M2 



2n 



M+Mm 



MMi~f 



M \ 1 



Mi^^P[-Wl]^^l^\<^\ (l-i^) 



M+Mn I 



Similarly, the bound-state contribution becomes 

I 1 



1 «:o(VM2 - + 



a;' 



/2 



^ - kI^MI, - 



Ml 



exp 



Ml u'-M 



(a^O), 
(a — >• 00), 



2^' + Mjv 

27r — c, 



where ur' 



M2 



Ma7 In the limit a — > cxo, 



r /I 
«o^— 1- — + - 
a L 2a 

a/ Va^ 



and the bound-state contribution is simplified as 



exp 



Ml, 



2.^ + ^2a(l 



2a/ Vay 



00 



MMjv 

Let us suppose that one can freely change the interaction strength of hyperon-nucleon 
and examine how the integral / should change as a function the interaction strength. When 
the interaction is weak, the scattering length is also small, and the integral / is dominated 
by/*: 



I = h + Ic 
1 



Ml 



exp 



M2 
Ml, 



27r- 



M + M, 



iV 



a + 0{a^). 



MMn 

As the interaction becomes stronger, the scattering length increases and the integral, /, also 
increases. As the interaction strength increases further, the scattering length eventually 
diverges when the bound state is just formed. Just before the bound state is formed, the 
integral / becomes 



I = h + Ic 

1 



Ml 



'"prMr"^ifM;r2+^U'' 



and just after the bound state is formed, it becomes 



I^It + Ic + h 



This shows that before and after the bound state is formed the integral is continuous, though 
the scattering length diverges with opposite signs. This observation leads us to conjecture 
that the integral around the threshold is a measure of the YN interaction strength. 
Here we study the integral in a solvable case, namely the separable potential model: 

We obtain 

1 / a 1 \ 



Ml ^ \ Ml) M/34 ' \M},)' 
We find that the integral / is proportional to the potential strength a if the range parameter 
(3 is fixed. 

Based on this conjecture we define the YN interaction strength, v, by 

1 ( M^\ 

In the dispersion integral, the imaginary part of the correlation function, Imn+, contains 
the contribution from all possible intermediate states such as those of the YN , YNtt channels 
and so on. However, only the YN channel contributes around the threshold. We assume 
that the contribution from the YN state is taken into account by the form of the right-hand 
side of Eq. (10) and that the rest is approximated by an asymptotic form of the correlation 
function starting from an effective threshold, uj+, for the B — 2 channels other than the YN 
channel and a;_ for the B — channels: 

(qo - M)ImH+ = X'^n5{qo - M)v 

+ {ei-qo - + e{qo - 0;+)} (qo - M)ImH^^^, (11) 

where H^^^ is the asymptotic form of the correlation function in the OPE and the normal- 
ization constant A is explicitly included {{0\'r]{0)\qs) — Xu{qs)). This is possible now because 
the contribution from states other than those of the YN channel is exponentially suppressed 
compared to the YN contribution. 
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III. OPE OF CORRELATION FUNCTIONS 



Let us turn to the OPE. We take the interpolating fields of hyperons [8] as 

f2 



for the A and 



for the other hyperons where qi — u, q2 — s for the E"*", qi — d, q2 — s for the E~, qi — s, 
q2 = u for the S° and qi — s, q2 — d for the S~. Since the interction strength for the 
I]^N channel is obtained by averaging those for the H'^N and T,~N, rj for the is not 
needed here. C denotes the charge conjugation operator and a, b and c are color indices. 
We take into account all the operators of dimension less than or equal to four. We also 
include four-quark operators of dimension six, since four-quark operators are known to give 
the largest contribution among higher order operators in the QCD sum rule for the nucleon 
mass [2,14]. The OPE of the correlation functions are given as 

+qnq''H-qV{ - ^(di^d) - ^{uiuu) - '^{s^us)} 

+g^ln(-g2)7r2{i ^^G^'^G,^^ + l(3md + - 2m,)(rfd) 
+^(3m„ + md- 2ms){uu) ^(Sm^ - 2md - 2m„)(ss)| 



+5'^ln(-gV{-^(-'5[Gj:G 



+ ^i{S[d^^D,d]) + ^i{S[u^^DM) + y^(5[s7M^.s])} 



+qt,-^'K^{^{ddss) ^{uuss) - ^{dduu)^ 

^[gMn(-g2){l(2m„ + 2ma - m,)(I)} + qHn{-q^)7:^[ - ^(Jd) - '^{uu) + ^(55)} 

1 - 1 

-Fg''ln(-5^)7r^|-(4md - 2m„ ms){d'y^d) + -(4m„ - 2m<i ms){u'y^u) 

2 _ ^ 
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116 16 
+q''—TT'^{ — {{ddu-ff,u) + (dds-fi^s)) + —{{uud-i^,d) + {uus^i^s)) 

q ^ y y 
--{{ssd^i_,d) + (55-^7^^))} 

+—n^^-mu{—2{uudd) + 12{ddss) + (uuss)) 
q y 

4 _ - 

+ -'md{—2{uudd) + 12{uuss) + (ddss)) 
9 

8 - 
+-ms{—{ddss) + 6{uudd) — (uuss)) 
9 

+^{2{uu) + 2{dd) - {ss}) (^^G^^G^fs) }' 
+ ^7''75[g^ln(-g2)7r2{^(J7^75d) + l-{u^^^^u) + y (57^755)} 



+q^q''ln{-q^)TT^{ - -{dj^-f^d) - -(^7^75^) - ^{siul5s)} 

(A - 4 20 

+g''ln(-g2)7r^|-(5[(i7^75iL>^(i]) + -(5['U7^75iL>^'u]) - — (5[s7^75iL>^s]) 



^(m<i + m„ - 2ms) {{djb'icrni.d) + {u-f^ia^,yu) - 2(575^(7^^5)) | 

116 16 

7r*^| - —{{ssd'^ir^ia^^d) + {ssu'^ir^ia^^u)) - —{{dds'^ir^ia^^s) + {uus'^ir^ia^ys)) 
q ^ \) y 

8 - - 
+-({ddu'y5ia^^u) + {uud'y^ia ^^d))^ 

1 111 

+ ^^75C^''''[g^ln(-g2)7r^{-((i75icT/,^ci) + - (M75icT^i.M) - ^(s75^o-^i.s)} 

+g^g''ln(-g^)7r^{^((i75i(7^pci) + ^ ('U75io-,.pii) - ^(s75«o-i.ps)} 
+g^ln(-gf^)7r^|^(3ms - Gm^ - Ama) {djul^d) + ^(Sm^ - 6md - ArUu) {u'yu^5u) 

--(3mrf + 3m„ - (57^755)} 

116 16 
+?/i^7r^lir((f^c?^^7t^75^i) + (dds-f^-f^s)) + —({uud-f^-f5d) + (uus-f^-f^s)) 
q^^ ^ y 9 
g 

--{{ssd'^y'^^d) + {ssu'^^'^r^u 

))}] (12) 

for A and 

n^^^(,) 

= i^7''[g.g'ln(-g^)33^(I) + q'H-<fV{ - ^(?-i7.gi) - ^(g27Mg2)} 

2 2 
+g^g''ln(-52)7r2{ - -(?i7i/?i) - ^{q2luq2)] 

+qM-<fy[^,Mq2) + ^ {^G-^G^p) } 

-rq^H-q'y{ - \ (^^S[G^^G + ^-^i{S[q,^,DM) + ^i(5[?-27^i?.?2]) } 
1 8 
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+ 



47r4 



+ 



ln(-g2) {Im,, ( I ) } + H-q'V{-{<l2q2)] 

1 16 

+g''ln(-g2)7r^{2mg2 ((527^52) - (gi7M?i))} + ^^^^^^{y (?2g2gi7M?i)} 
+ ^7r^{8m<,i(gigig2g2) + ^rnq^{Mm(li) + ^(^292) {^G'^'^Ga^ } 

4^7^75 [g^ln(-g2)7r^{^(gi7^75gi) - ^ (g27/x75g2) } 

2 2 
+g^g''ln(-g^)7r^{ - -(gi7i/75gi) - ^ (?27i/75?2) } 

8 4 
+g''ln(-g2)7r2{ - -{S[qi'^^'^r,iD^qi]) + -(5[g27/.75^^i.g2])} 

+g''^7r^{ - y(gigigi75«o-/.i.?i>} 

11 2 

+ 4^^75C^^''[g^ln(-g^)7r2{-(g275^cr^i/g2)} + q^lq''^'^{-q^)^'^{^{q2lb■i(Tupq2)}} 

+g^ln(-g^)7r^{ - -mg2((g27i^75g2) + 3(gi7^75gi))} 
+?M4^^{y(^2g2gi7i/75?i)} 



(13) 



for other hyperons where the matrix elements can be taken with any state. The vacuum-to- 
vacuum and nucleon-to-nucleon correlation functions, n^'^^(g) and n'^'^^(gps), are obtained 
by replacing the matrix element of (O), with the vacuum and nucleon matrix elements, 
(0)0 = (0|O|0) and {0)n = {ps\0\ps) - {ps\ps){0)o, respectively. 

We now discuss the nucleon matrix elements of the quark-gluon operators. Dimension- 
three operators are quark bilinear operators, qq, q^^q, Tf^'^^q and q^Qcr^uq- The nucleon 
matrix elements, (qq) n and (q'Jij.q) n, are spin-independent and have already been discussed 
in Ref. [9], while {q'~i^/)^q)N and {q'~i50'^yq) jq are spin-dependent and are written in terms of 
the axial charge, Ag, and the tensor charge, 5q, as 

{qinibq)N = Ags^, 

{qilb'7iivq)N = 5q{s^p^ - s^p^)/po- 

For the vector and scalar charges, we take, {u^u)p — 2, {cfd)p — 1, {s^s)p — 0, {uu)p — 3.46, 
{dd)p = 2.96, {ss)p = 0.7 [10]. For the axial and the tensor charges, we use the values, 
Au = 0.638 ± 0.054, Ad = -0.347 ± 0.046, As = -0.109 ± 0.030, Su = 0.839 ± 0.060, 
5d = -0.231 ± 0.055, Ss = -0.046 ± 0.034, obtained by the lattice calculations [11,12]. 
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Dimension-four operators are gluon operators, ^G'^'^G^j,^ ^^SfG^Gp,^], and quark oper- 
ators, S[q'^^iDyq\, qS{^i^iDi,Yi^q. The nucleon matrix elements of the gluon operators are 
spin-independent and have already been discussed in Ref. [9,13] 

(^G^^G^p = -738 MeV, (^5[G'^oG"^°])p = -50 MeV. 

TT TT 

The matrix elements of the spin- independent operators, (S[g7^iZ^i,g], have also been discussed 
in Ref. [14] 

i{S\u-i^D^u])p = 222 MeV, i{S[d-fi,D^d])p = 95 MeV, 
i{S[s-f^D^s])p = 18 MeV. 

On the other hand, the matrix elements of the operator, qS{'yfj,iDi,)j^q, is spin-dependent 
and is expressed as 

{qS{-f,^iD^)-f5q)N = aiS{^p,.}, 

where ai is related to the first moment of the longitudinal quark-spin distribution, g2, if 
operators including the quark mass are neglected, which is given at the tree-level by [15] 

ai — —2 / dxxg2(x). 
Jo 

Using this expression we have calculated oi from g2. Experimental data for g2 are given in 
Ref. [16] over the range 0.075 < x < 0.8 and 1.3 < < 10 ((GeV/c)^) for the proton and 
in Ref. [17] over the range 0.06 <x < 0.70 and 1.0 < < 17.0 ((GeV/c)^) for the neutron. 
The results are as follows: 

< = 0.05 ±0.04, of = -0.08 ±0.13. 

Unfortunately a* is not determined by experimental data up to now. However, absolute 
values of nucleon matrix elements of strange quark operators are generally much smaller 
than those of up and down quark operators. Therefore, we set af = 0. 

The dimension-six four-quark operators are qqqq, qqTi^iq-, qqqi/xIbQ and qqq'J^cr^uQ- The 
nucleon matrix elements of the four-quark operators are approximated to factorize [7] as 

14 



{qQQQ)n = 2(gg)o(gg)iv, {qqq^nq)N = {qq)o{qinq)N, 

{qqqinim)N = (^g)o(^7M75g)iv, {qqqib(7^ivq)N = {qq)o{qib(^t,vq)N, 



which can be evaluated from the vacuum and nucleon matrix elements of dimension-three 

quark operators. 

For completeness we also list here the values of vacuum condensates and quark masses 

{uu)q = {dd)o = -(250 MeV)^ {ss)o = -(232 MeV)^ 
{^G^)o = (330 MeV)^ 

TT 

mu^irid^ 0, rris = 120 MeV. 



IV. RESULTS 

The OPE of the projected correlation function, 11^^^, for the A or other hyperons is given 
by substituting Eq. (12) and Eq. (13), respectively into the right-hand side of Eq. (3). Then, 
substituting 11^^^ into the left-hand side of Eq. (7), and Eq. (11) into the the right-hand 
side of Eq. (7), respectively, we obtain the Borcl sum rules for hyperon-nucleon interaction 
strengths. The sum rule for the spin-independent part is 



exp 



M2 
Ml, 



y2^indep. 



47r4 



{ (CsMMb - C3MI) 

{ 



D o 3 

+ {C^M - C^Mb) [^{\ l^G'^^G. 



'■{dd) 



N 



- -fa/3/ - \—<S[GoGpq] 

8 \ TT In oXtt in 

16 - 16 28 

+—i{S[d'-foDQd])N + —i{S[u'-fQDQu]) N + —i{S[s'-iQDQs\) ^ 

-F^(3md m„ - 2ms){dd)isf + ^(3m„ + md- 2ms){uu)N 

-F^(3ms - 2m<i - 2mu){ss)N + ^(4md - 2m„ + ■ms){djod)N 



1 2 



M 



Ml 



1 6 

TT^I —{ddss)N + -^{uuss)n — -z{dduu)N 
^9 9 9 



16 



16 - _ _ 16 _ - 

+ — {dd)Q{{ujou)N + (s7o-s)iv) + y (MM)o((c?7oc?)Ar + {sJqs)n) 



15 



g 



(14) 



for the AA^" channel and 



exp 



M2 



N 



— { {C2MMB - C^Ml) [tt^I - 3(glgi)jv - (gk2)iv} - ^^^^2) 

'-'^OXTT /AT 0\7r /AT 

16 _ 4 _ 

+ yi(cS[gi7oDogi])iv + 2«(5[g27o^og2])jv 

16 



+ 



M 



Ml 



^7r''(gigi?igi)iv + ^7r''(g2g2gigi)iv] } 



(15) 



-3 " 3 

for other hyperon-nucleon channels. Similarly, the sum rule for the spin-dependent part is 



Ml 



exp 



Ml, 



y2^dep. 
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{ss)o{{d-fkl5d)N + (M7fc75M)iv)} 



(16) 



for the AA" channel and 
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y7r^(?i?i?i75«o-fcogi)iv - y7r^(?2g2?i7fe75?i)iv] } 
for other hyperon-nucleon channels. Ci, C2 and C3 are given by 



(17) 



Ci = l 



exp 



Ml, 



+ exp 



Ml, 



exp 



M, 



Ml, 



C. = l- 



Mb 
2 



exp 



+ 



1 + TTT I exp 



^1. 



^1, 



1 



L0_ 

Ml. 



exp 



a; 



Mb 
2 \ " 



- $ 



^1. 



where $ is the error function. 

Prom the vacuum-to- vacuum correlation function we obtain the Borel sum rule for the 
normalization constant as [2,8] 
1 ( 
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exp 
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47r4 



1 1 
+ -(3m„ + md- 2ms){uu) + -(Sm^ - 2m^ - 2m„)(ss)| 

1 4fl6,7,_, 16 , _ , 8,y, ,v 



for the A and 
1 
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exp 
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1 



47r4 



- D,Mi\ - (^G'-^'G,;,) +m,,(g2?2)o} - -^{^'\{Mimi)o\\ (19) 



M| I- 3 



for other hyperons. D\ and D2 are expressed by the effective threshold, ujq, as 



(18) 



Di = 1 — exp 



i^, = i-(i+ ' 



Ml 2Ml^ 



exp 



Dividing Eq. (14) (Eq. (16)) and Eq. (15) (Eq. (17)) by Eq. (18) and Eq. (19), respectively, we 
finally obtain the sum rules for the spin-independent (spin-dependent) part of the hyperon- 
nucleon interaction strengths. 
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It is instructive to look at sum rules in the leading order of the OPE before discussing 
the full results: 
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{5.46 + s-s' [0.674]}, 

1 1 
{3(s^s)p + {dU)p + {dd)p + s • s ' [-(-5As + Ad) + -5(i] } 

{3.96 + s-s' [-0.011]}, 



where s (s ') denotes the hyperon (nucleon) spin operator. 

For comparison, the interaction strength in the neutron-proton {np) channel is given by 

{?>{d^d)p + {v)u)p + {uu)p + s • s' ^(-5 Ad + Au) + ]-5u } 



~ Ml 

Stt^ 



{8.46 



s • s 



1.07 



}• 



The spin-independent interaction strength is related to the vector and scalar charges of 
the nucleon while the spin-dependent strength to the axial and tensor charges. In the A^A^ 
channel only up and down quark operators contribute, while in the YN channels strange 
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quark operators appear in addition. It is interesting that the ratio of different flavor con- 
tributions is the same for the scalar and tensor charges in all YN and A'"A^ channels. If 
the proton matrix elements of the up, down and strange quark operators are the same, all 
YN interaction strengths coincide, which holds even if higher order terms of the OPE are 
included. It should be noted, however, that this is not the SU{3) hmit. Even in the SU(3) 
limit the proton matrix elements of quark operators with different flavors are different in 
general. 

Inserting the numerical values for the matrix elements we find that the spin-dependent 
interaction strength is considerably smaller than the spin-dependent interaction strength 
in each channel. This is because the absolute values of the axial and tensor charges are 
considerably smaller than the vector and scalar charges except for the strange vector charge, 
which identically vanishes. Also found is that both the spin-independent and spin-dependent 
interaction strengths are smaller in the YN channel than in the NN channel, since the 
strange charges are much smaller than the up and down charges. An exception for this is 
the spin-dependent interaction strength in the E+p channel, which receives large contribution 
from the up axial charge. Among YN channels the AN has the largest interaction strength 
due to the largest scalar charge contribution, which is even larger than in the NN channel 
since the sign of the Wilson coefficients is oposite to the others. If one looks at the charge 
dependence, the interaction strength is greater in the E+A?" (EPN) channel than in the E~A^ 
(S~A^), since the up charges are larger than the down charges. 

Let us now discuss the results including all the operators up to dimension 4 and the 
dimension 6 four-quark operators. Up to now we have not explained how to determine cuq, 
UJ+ and CU-. They are determined by the Borel stability analysis, i.e. , so that the calculated 
interaction strength has the most stable plateau as a function of the Borel mass squared. We 
constrain a;_ to be greater than cuq since we expect the projection reduces the continuum 
contribution in the negative low energy region. Also, we take a;+ to be equal to cuq for 
simplicity. The interaction strength, however, is not sensitive to the choice of cu^. Even if 
we take |a;o — u;+\ — 0.5 GeV, the interaction strength changes only less than 3%. Fig. 1 
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shows the calculated interaction strength as a function of the Borel mass. The upper and 
lower curves beyond = 1.2 GeV^ correspond to the spin-independent and spin-dependent 
interaction strengths, respectively. There appears a stale plateau above M| = 1.4 GeV^ in 
all channels. The prediction of the Borel sum rule is determined by taking the value in 
the region of the plateau. The optimal value of the continuum threshold, cu+, in the E~p 
channel is much different from those in other channels. However, even if we take a;+ in the 
E"^p channel to be the same as in the E^p channel, the difference in the interaction strength 
is less than 10% for M| = 1.0 ~ 2.0 GeV^. Fig. 2 shows the calculated interaction strengths 
without including the continuum. From Fig. 1 and Fig. 2 we find that the contribution of 
the continuum is rather small. 

Comparing the results in the leading order of the OPE and those including the higher 
order terms we find that the contribution of the higher order terms makes the strength 
decrease in the UN channel but increase in the EN channel, relative to the strength in the 
AN channel, so that the E+p strength becomes smaller than the E^p strength. We see that 
the spin-independent interaction strength in the NN channel is greater than those in the 
YN channels, which is similar to the results in the leading order. We should, however, keep 
in mind uncertainties in the matrix elements and the continuum contributions and therefore 
conclude at the moment that the strengths in the np, Ap and E^p channels are similar. The 
strengths are smaller in the Ep and E^p channels than in other three channels. Furthermore, 
the E~p interaction is the weakest in all the hyperon-proton interactions. 

Comparing the spin-independent part with the spin-dependent part we find that the 
spin-dependent interaction is much smaller than the spin-independent interaction in each 
channel. The interaction strengths in the E+p and E~p channels have an opposite sign to 
the interaction strengths in other channels. The spin-singlet interaction is greater than the 
spin-triplet one in the E"'"p and E~p channels but the relation is opposite in other channels. 

Up to now we have discussed the results of sum rules identifying them as interaction 
strengths because interaction strengths are convenient for comparing different hyperon- 
nucleon and nucleon-nucleon channels. Now, we want to compare our results with those 
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of phenomenological meson-exchange models. For this purpose, scattering lengths are more 
convenient. Therefore, from now on we identify our results as scattering lengths ignoring 
0{a^) differences. Table I shows calculated scattering lengths with estimated errors for 
the hyperon-nucleon channels in fm. The errors due to uncertainties of the spin-dependent 
nucleon matrix elements are as small as ±0.04 ~ ±0.07 fm. Among the spin-independent 
nucleon matrix elements the scalar charges cause most of errors. The uncertainties of {uu)p 
and {dd)p are about 30% and that of {ss)p is about 60%. These errors change the scattering 
lengths by ±0.4 ~ ±1.2 fm. {S[q^ij,D„q])p is expected to have much smaller error than 
{qq)p [18]. The uncertainty of {^G ^1,0^'') p is about 20%, which changes asp by only about 
2% and affects others less than 0.3%. The errors of other spin-independent nucleon matrix 
elements are also thought to be small since they are higher order terms. Furthermore, the 
vacuum condensates of the quark and the gluon operators have uncertainties about 20% 
and the strange quark mass lies in the range between 100 MeV and 200 MeV. The total 
error due to these uncertainties is about 15%. The errors of the scattering lenghts due to 
the uncertainties of all these input parameters are shown in Table I. Even though we did 
not include in the errors of Table I, we also estimated errors due to effective continuum 
thresholds which change oy^r roughly by ±1 fm. 

Table 11 shows the scattering lengths calculated in the Nijmegen meson-exchange models, 
NSC97a-f [19]. Comparing Table I and Table II, we find that the scattering lengths obtained 
by sum rules are rather larger than those in the meson-exchange models. It should be 
remembered, however, that the results of sum rules coincide with the scattering lengths only 
when the scattering lengths are small. The results of the sum rules and the meson-exchange 
models have some common tendencies: the E+p channel has stronger spin dependence than 
the Kp channel and the triplet scattering length is smaller than the singlet scattering length 
in the E+p channel. On the contrary, the results are quit different from each other in the 
S°p channel. This might be related to anormalous behavior of the effective range in the 
Nijmegen models, which is extremely large in the triplet channel and is negative and large 
in the singlet channel. 
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A comment on the spin dependence of the interaction in the AA^ channel and observed 
hypernuclear states is in order here. Experimentally, the ground state of ^He-hypernucleus 
is known to be 0"*" and 1"*" is an exited state. This seems to contradict with the present 
result that the triplet Ap interaction is more attractive than the singlet one. According to 
Akaishi et al., however, if one takes the three-body force into account the ordering of the 
observed ^He-hypernuclear states is reproduced with the Nijmegen model D [20] in which 
the triplet scattering length is a httle greater than the singlet one [21]. Even though this is 
not conclusive, this at least tells us that one must be careful when one extracts information 
on the hyperon-nucleon interaction from observed hypernuclear states. 



V. SUMMARY 

In summary, we have formulated the QCD sum rule for the YN interactions starting 
from the spin-dependent hyperon correlation function, whose matrix element is taken with 
respect to the one-nucleon state. The dispersion integral around the threshold is identified 
as a measure of the interaction strength in the YN channel. The Wilson coefficients of the 
OPE for the hyperon correlation functions have been calculated. Strange quark operators 
appear as well as up and down quark operators. After Borel transforming we have obtained 
Borel sum rules which relate the YN interaction strengths to the nucleon matrix elements 
of quark-gluon operators. We have then discussed physical implications of the obtained sum 
rules for the interaction strengths in the YN and the A^A^ channels. It was found that the 
interaction strength is smaller in the YN channels than in the A"A^ channel since the matrix 
elements of strange quark operators are smaller than those of up and down quark operators. 
Among YN channels the AA" channel turned out to have a stronger interaction strength 
than the EA" and SA" channels. It was also found that the E+A" (EPN) interaction strength 
is greater than the E~A^ ['E~N) interaction strength since matrix elements of up quark 
operators are greater than those of the down quark operators. Taking uncertainties into 
account, however, we tentatively concluded that the np, Ap and S°p channels have similar 
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interaction strengths. It turned out that the strengths are smaller in the Ep and S~p 
channels than in other channels and also that the channel has the weakest interaction 
in all channels. Also found is that the spin-dependent strength is smaller than the spin- 
independent strength in all channels. The spin-dependence in the E+p and E~p channels 
is oposite to other channels. The results are compared with those of the phenomenological 
meson-exchange models. The results of the sum rules and the meson-exchange models have 
some common tendencies. They are, however, different in some other respects. 

In the present paper we have ignored the coupling of different hyperon-nucleon and 
hyperon-hyperon channels such as AN-T^N and SA'^-EE-AA. One might be able to formulate 
the sum rule for the transition, AN — > EA^, by considering the time-ordered product of the 
sigma field and the lambda field, which is an interesting extension of the present approach. 
Another interesting possibility is to study hyperon matrix elements of the quark-gluon com- 
posite operators by means of the sum rule. This would be done by constructing sum rule 
starting from the hyperon matrix element of the nucleon correlation function and relat- 
ing the obtained hyperon-nucleon interaction strengths with the ones in the present paper. 
These applications would provide us with important new information about hyperon-nucleon 
interactions and hyperon structures. 
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TABLES 

TABLE I. The calculated scattering lengths for the hyperon-nucleon channels (uint in fm). 
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I 

> 

O 



m 
> 



400 



200 



-200 



-400 




J I I L 



J I I L 



J I I L 



J I I L 



0.5 



1.0 



1.5 



2.0 



2.5 



U/ (GeV^) 



FIG. 1. The interaction strength, VBp, as a function of the Borel mass where ujq = 2.8 GeV 
and u- = 3.4 GeV for the np channel, = 2.7 GeV and a;_ = 3.3 GeV for the Ap channel, 
a;o = 3.1 GeV and a;_ = 4.0 GeV for the S+p channel, = 3.5 GeV and w_ = 10.8 GeV for 
the channel, loq = 2.4 GeV and lo- = 2.8 GeV for the S~p channel and ujq = 2.5 GeV and 

= 2.9 GeV for the E~p channel. The np, Ap, S+p, T>~p, S+p and E~p channels are shown as 
solid, dotted, dashed, thick dashed, dot-dashed and thick dot-dashed lines, respectively. 
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FIG. 2. The interaction strength, VBp, as a function of the Borel mass where (Jq = oo and 
u± = oo. The np, Ap, S+p, T,~p, S+p and 'E~p channels are shown as sohd, dotted, dashed, thick 
dashed, dot-dashed and thick dot-dashed hnes, respectively. 
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